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Abstract 

We derive explicit forms of Markovian transition probability densities for the velocity space, 
phase-space and the Smoluchowski configuration-space Brownian motion of a charged particle in a 
constant magnetic field. By invoking a hydrodynamical formalism for those stochastic processes, 
we quantify a continual (net on the local average) heat transfer from the thermostat to diffusing 
particles. 

1 Introduction 



We address an old-fashioned problem of the Brownian motion of a charged particle in a constant 
magnetic field. That issue has originated from studies of the diffusion of plasma across a magnetic 
field j^, 1^ and nowadays, together with a free Brownian motion example, stands for a textbook 
illustration of how transport and auto-correlation functions should be computed in generic situations 
governed by the Langevin equation (to a suitable degree of approximation of a kinetic theory, when 
collisions are stochastically modeled in terms of a random force), cf. but also Q, 

From a purely pragmatic point of view this white-noise strategy is quite satisfactory. After (for- 
mally) evaluating velocity auto-correlation functions, formulas for running and asymptotic diifusion 
coefficients easily follow. To that end an explicit form of the probability density or transition probabil- 
ity density of the involved stochastic diffusion processes (in velocity space, phase-space or configuration 
space) is not necessary, cf. [Q, 

To our knowledge, except for the paper (mentioned in ||] as a footnote reference for the purpose 
of evaluation of the mean square velocity and its mean square displacement at equilibrium), for a 
Brownian particle in a constant magnetic field no attempt was made in the literature to give a complete 
characterization of the stochastic process itself, nor pass to the associated macrosopic (hydrodynamical 
formalism) balance equations. (Cf. ||], 0, [||, [0, for a number of reasons why to do that). 

Surprisingly enough, in Ref. the Brownian motion in a magnetic field is described in terms of 
operator-valued (matrix-valued functions) probability distributions that additionally involve fractional 



1 



powers of matrices. In consequence, there is no clean path towards a (necessary) relationship with the 
associated Kramers- Smoluchowski equations (cf. Chap. 6.1 in Ref. Q), nor ways to stay in conformity 
with the standard wisdom about probabilistic procedures valid in case of the free Brownian motion 
(Ornstcin-Uhlenbcck process), cf. |]. 

Therefore, we decided to address an issue of the Brownian motion in a magnetic field anew, to 
unravel its features of a fully-fledged stochastic diffusion process. In particular, we derive transition 
probability densities governing both the velocity, phase-space and the configuration space processes. 
Hydrodynamical balance equations and their behaviour in the Smoluchowski regime are discussed as 
well. 

2 Velocity-space diffusion process 

The standard analysis of the Brownian motion of a free particle employs the Langevin equation 
= —pijt+ A (t) where denotes the velocity of the particle and the influence of the surrounding 
medium on the motion (random acceleration) of the particle is modeled by means of two independent 
contributions. A systematic part — /3l? represents a dynamical friction. The remaining fluctuating 
part A {t) is supposed to display a statistics of the familiar white noise: (i) A (t) is independent of 
~Tt, (ii) {Ai (s)) — and {Ai (s) Aj (s')) = 2qSij6 (s — s') for i,j — 1, 2, 3, where q = is a physical 

parameter. The well-known Ornstcin-Uhlenbeck stochastic process comes out on that conceptual 
basis, |l|,ll|, g. 

The linear friction model can be adopted to the case of diffusion of charged particles in the presence 
of a constant magnetic field which acts upon particles via the Lorentz force. The Langevin equation 
for that motion reads: 



^ = -/?lf+^lI>x 5^ + ^(0 (1) 
at mc 

where qe denotes an electric charge of the particle of mass m. 

Let us assume for simplicity that the constant magnetic field is directed along the z-axis of a 
Cartesian reference frame: B = (0, 0, B) and B — const. In this case Eq. (|^) takes the form 

— ^-AU^+A{t) (2) 

where 

f P -iOc \ 
A = \ ij, P \ (3) 
V p J 

and iUc = denotes the Larmor frequency. Assuming the Langevin equation to be (at least 
formally) solvable, we can infer a probability density P {Ut, t\lto), t > conditioned by the the initial 
velocity data choice it = Ttg at t = 0. Physical circumstances of the problem enforce a demand: (i) 

P{lj^,t\Tta) 6^{lt ~ Tta) as i -> and (ii) P{lt,t\lto) [ 27,Tbt ) ' ''^P (" "LbT ) as t ^ oo. 
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A formal solution of Eq. (m) reads: 



Tt (t) - e-^'Tto = I e-^^'-''> A (s) ds . (4) 
By taking into account that 







e 



cosojc^ sinwc^ 
-At^g-/3tj _gintj^t cos Wet \=e''^*U{t) (5) 
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we can rewrite (||) as follows 



(t) - e-^*C/ (<) H^o = / e-^^*"'^) U{t-s)A (s) ds . (6) 





Statistical properties of (t) — e~^*ir'o are identical with those of A {s)ds. In consequence, 
the problem of deducing a probability density P(H^,i|ir'o) is equivalent to deriving the probability 
distribution of the random vector 

5"= / V(s)^(s)ds (7) 
Jo 

where i/' (s) = e-^^*-") e-^(*-^)C/ (i - s). 

The white noise term A (s) in view of the integration with respect to time is amenable to a more 
rigorous analysis that invokes the Wiener process increments and their statistics, |l4|. Let us divide 
the time integration interval into a large number of small subintervals At. We adjust them suitably 
to assure that effectively (s) is constant on each subinterval {jAt, {j + 1) At) and equal ip (jAt). As 
a result we obtain the expression 



N-l Aj+i)^t 
^=Y.iP{jAt) A{s)ds. (8) 

J]At 



3=0 



Here (At) = Jj j^^^'^^ ~^ (s) ds stands for the above-mentioned Wiener process increment. Phys- 
ically, B (At) represents the net acceleration which a Brownian particle may suffer (in fact accumu- 
lates) during an interval of time At. 

Equation (||) becomes 

Af-l JV-1 

^ =J2^{jAt)V{At)^Y.^^ (9) 

i=o j=o 

where we introduce 's^j — ■0 (jAt) (At) ~ i^jl^ (At). 
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The Wiener process argument [|l2[ , , g allows us to infer the probability distribution of ~s^j . It 
is enough to employ the fact that the distribution of (Ai) is Gaussian with mean zero and variance 
q = Then 



B (At) 



( 1 



exp 



(10) 



and in view of j — ^pj 1? {At) by performing the change of variables in (|T^) we get 



=det [^p-^]w[ijT^^^] = 



det -0,- 



Since det V' (s) = e-3/3(t-^) and (s) ^ U [- (t ~ s)] e^^*"^' we obtain 



w[s j\ 



AnqAt 



1 



-3/3(t-jAt) 



exp 



|g/3(t-jAt)^[_ (i-jAt)] 

4qAt 



and finally 



/ 1 



1 



exp 



4gAte-2/5(t-jAt) 



(11) 



(12) 



(13) 



Clearly, "s*^- are mutually independent random variables whose distribution is Gaussian with mean 
zero and variance tT| = 2qAte~'^^^*~^^^\ Hence, the probability distribution of = SjLo^ ~^ j 
again Gaussian with mean zero. Its variance equals the sum of variances of j i.e. = ct| = 
2qY.^Ate-^P^'-3^'\ 

Taking the limit N ~^ oo {At 0) we arrive at 



Jq 



(1 - e~2'3*) 



(14) 



Because of if the transition probability density of the Brownian particle 

velocity, conditioned by the initial data ll^o at <o = reads 



,-At 



P{'Tt,t\TtQ) 



271 



(1 



-2()t 



exp 



u 



2hBT ,^ _ g_ 



(15) 



The process is Markovian and time-homogeneous, hence the above formula can be trivially ex- 
tended to encompass the case of arbitrary t^ ^ Q : P (H*, tjir'oj ^o) arises by substituting everywhere 
t — to instead of t. 
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Physical arguments (cf. demand (ii) preceding Eq. (4)) refer to an asymptotic probability distri- 
bution (invariant measure density) P{u) of the random variable ~it in the Maxwell-Boltzmann form 

"•^''IsS^rj «»(-^j ■ <"» 

This time-independent probability density together with the time-homogeneous transition density 
(15) uniquely determine a stationary Markovian stochastic process for which we can evaluate various 
mean values. 

Expectation values of velocity components vanish: {ui {t)) = UiP {it) dlt = for i — 1,2,3. 
The matrix of the second moments (velocity auto-correlation functions) reads 



{ui (t) Uj {to)) = UiU^P {lf,t;Tto,to) dltdlto (17) 

J —oo 

where i,j = 1,2,3 and in view of P (IT", H^o, ^o) = ^ (^, ^|^o, ^o) ^' (^o) we arrive at the 
compact expression 



iy\t-ta\ _ ^bT_ l3\t-to 



costJc 1^ — ^o| sinwcl^ — ^o| 



g-A|t-to| ^ ^_g-P|t-to| -smu;c\t-to\ cosLJc\t-to\ . (18) 

\ Q 1 

In particular, the auto-correlation function (second moment) of the x-component of velocity equals 

(ui {t) ui {to)) = M^e-'^l*-*"! coso;, \t ~ to\ (19) 

TO 

in agreement with white noise calculations of Refs. and cf. Chap. 11, formula (11.25). 
In particular, the so-called running diffusion coefficient arises via straightforward integration of the 
function i?ii(r) =< ui{t)ui{to) > where t — t ~ to > 0: 

* ksT P + [LLicSin{ijJct) — (3cos{uJct)]exp{— pt) 

TO P'^ + UJc 

ksT 



with an obvious asymptotics (the same for D2{t)): Db ~ limt^ocDi{t) — ^T" ^ large 

fcflT 
mf3 



friction {lJc fixed and bounded) version £) — 



3 Spatial process 

The cylindrical symmetry of the problem allows us to consider separately processes running on the 
XY plane and along the Z-axis (where the free Brownian motion takes place). We shall confine 
further attention to the two-dimensional Xy-plane problem. Henceforth, each vector will carry two 
components which correspond to the x and y coordinates respectively. We will directly refer to the 
vector and matrix quantities introduced in the previous section, but while keeping the same notation, 
we shall simply disregard their z-coordinate contributions. 
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We define the spatial displacement of the Brownian particle as folows 



T* — T*o = / iv) dn (21) 

where ~li (t) is given by Eq. (|^) (except for disregarding the third coordinate). 
Our aim is to derive the probability distribution of at time t provided that the particle position 
and velocity were equal ~r*o and U^o respectively, at time to = 0. 

To that end we shall mimic procedures of the previous section. In view of: 



/ e-^^M^o^ / drj ^ dse-^^"-'^^ (s) (22) 

we have 



Jo Jo 



- T^o - (1 - e-^*) ^0 = ^ (l - e^(^-*)) 'X (s) ds (23) 

where 



is the inverse of the matrix A (regarded as a rank two sub- matrix of that originally introduced in 
Eq. (3)). Let us define two auxiliary matrices 

n = A-i(l-e-^*) (25) 

0(s) ^ A-i (l - e^(^-*)) . 

Because of: 

-At _ P W _ -i3t ( cosujct sin Wei \ _ -pti 



e-^* = exp ^ T e-l'' -—c^ — ^ e-"'U {t) (26) 

[ \ ujc P J ) \ — smwc^ coswc^ 

we can represent matrices 17, (s) in more detailed form. We have: 



n ^ If ^ ) - e-^^ { ^ ^^"'"^'^ ) ^ (27) 

^2+tj2 |^\^ ^ y y P J \ -smt^ci cost^ci ' ' 

and 



0(s) = A-i(l-e-'3(*-^)f/(i-s)) = (28) 
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/3 iOc\ f 1 - e^(^-*)cosa;c(s-t) -e^^^^"*) sinwc (s - 



P'^ + oj'^ \ -ujc f3 J \ e^(«-*)sina;c(s-t) 1 - e'^^*"*) coswc (« - i) 

Next steps imitate procedures of the previous section. Thus, we seek for the probability distribution 
of the random (planar) vector ^ = Jq (s) ^ (s) ds where ^ = — T*o — filfo- 

Dividing the time interval (0, t) into small subintervals to assure that ^ (s) can be regarded constant 
over the time span {jAt, {j + 1) At) and equal (j) {jAt), we obtain 

Af-l Aj+l)At N-1 N-1 

where T'j = (jAi) ^ (At) = (j)jT^ [At). 

By invoking the probability distribution (10) we perform an appropriate change of variables: = 
{At) to yield a probability distribution of 

w [f^j] = det [0-1] w [</.- V,] = [ct>fT'j\ . (30) 

Presently (not to be confused with previous steps (11)-(15)) we have 

det cl> is) = ^„ ^ , f 1 + e''^^^-*) - 2e^^'-*^ cosue {s - t)) (31) 

and 



<^"^ (s) = ^7?^^^ 7 7 [l - e^(«-*)?7 (- (s - t))] A . (32) 

^ ^ l + e2/3(^-*) -2p/(^-*)coswc(s-i) L V V v ; 

So, the inverse of the matrix has the form: 



where 



(33) 



l_e/50-^*-*)cosu;e(jAt-t) e^^^*"*) sina;^ (jAt - f) \ / /? -w,. 
-e^(j^*-')smWc{jAt-t) 1 - e^^*"*) coswe (jAi - i) V '^c /? ' ^ ' 



and 



jj = l+ e2/30-A*-t) - 2e'5(i^*-*) cos We (jAt - t) . (35) 

There holds: 



7 



Ti 

and as a consequence the probability distribution of ~r^j becomes 



(36) 



w\r ■i\ = 



3^7j V47rgAt 



exp < 



(37) 



In view of 



that finally leads to 



A, 



(38) 



(39) 



Since this probability distribution is Gaussian with mean zero and variance (t| = 2qAt ^■^^^^ 'jj, 
the random vector ^ as a sum of independent random variables has the distribution 



w 



(7?) 



■ exp 



(40) 



In the limit of At — > we arrive at the integral 



(41) 



ds 



with Jq 7 (s) ds = t + 0, where 



(42) 



e = e{t) = ^{l-e-'f'')-2—^[P + {ucSinijj,t-PcosuJct)e-l''] . (43) 

That gives rise to an ultimate form of the transition probability density of the spatial displacement 
process: 
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with = r2(i) defined in Eqs. (|25|), (27). Notice that an asymptotic diffusion coefficient Db = 
■D ^a^^a of Section 3 (cf. Eq. (20)) appears here as a spatial dispersion - attenuation signature (when 
Uc grows up at /3 fixed). 

The spatial displacement process governed by the above transition probability density surely is not 
Markovian. That can be checked by inspection: the Chapman-Kolmogorov identity is not valid, like 
in the standard free Brownian motion example where the Ornstein-Uhlenbeck process induced (sole) 
spatial dynamics is non-Markovian as well. 

4 Phase-space process 

4.1 Free Brownian motion, Kramers equation and local conservation laws 

We take advantage of the cylindrical symmetry of our problem, and consider separately the (free) 
Brownian dynamics in the direction parallel to the magnetic field vector, e.g. along the Z-axis. 

That amounts to invoking a familiar Ornstein-Uhlenbeck process (in velocity /momentum) in its 
extended phase-space form. In the absence of external forces, the kinetic (Kramers-Fokker-Planck 
equation) reads: 



dtW + uV^W ^ f3\7uiWu) + qAuW (45) 

where q = D(3^. Here (3 is the friction coefficient, D will be identified later with the spatial diffusion 
constant, and (as before) we set D = ksT/mP in conformity with the Einstein fluctuation-dissipation 
identity. 

The joint probability distribution (in fact, density) W{z, u, t) for a freely moving Brownian particle 
which at i = initiates its motion at xq with an arbitrary inital velocity uq can be given in the form 
of the maximally symmetric displacement probability law: 

W{z,u,t) = W{R, S) = [An'iFG - H^T^I^ ■ e.p{_^^l^|^^±Z^} (46) 

where z- uq{1 - e^''*)/?-!, S = u- uoe"'^* while F = ^(2/3t - 3 -f 46"'^* - e'^''*) G = 
D(3{1 - e-2/3t) and H = D{1 - e-^^^f. 

For future reference, let us notice that marginal probablity densities, in the Smoluchowski regime 
(take for granted that time scales /3~^ and space scales (_D/3^^)^/^ are irrelevant [|2|) display familiar 
forms of the Maxwell-Boltzmann probability density w{u, t) — (2^^)^^^ exp{— ^^t ^ ^^"^ ^^"^ diffusion 
kernel w{z,t) = (47r£>t)~-^/^ea;p(— 2^) respectively. 

A direct evaluation of the first and second local moment of the phase-space probability density 
gives 
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<u>= I duuW{z,u,t) = w{R)[{H/F)R + uoe-^'] (47) 



<^'>^J duu'W{z,u,t) - ( ^^^^' + |Ji?2) . (2nF)-^'\xp{~^) . (48) 

Let us notice that after passing to the diffusion (Smoluchowski) regime, 0|, one readily recovers 
the local (configuration space conditioned) moment < u>z~ < u > to be in the form 

z ^\7w(z,t) 
<u>.^- = ~D / / 49 
2t w[z,t) 

while for the second local moment < > z= — < > we would arrive at 



<u^ >,= {DP-D/2t)+<u>l . (50) 

By inspection one verifies that the transport (Kramers) equation for W{z,u,t) implies local con- 
servation laws: 



dtW + \/{<u>zw) =0 (51) 

and 

dti< u>z w) + Vz(< >z w) = -P < u >z w . (52) 

At this point (we strictly follow the moment equations strategy of the traditional kinetic theory 
of gases and liquids, compare e.g. ||ll|) let us introduce the notion of the pressure function Pkm'- 

Pkiniz, t) = {<u^ <u >l)w{z, t) (53) 

in terms of which we can analyze the local momentum conservation law 

idt+ <u>zV)<u>,^~[3<u >z -Yfh^ . (54) 

w 

One should realize that in the Smoluchowski regime the friction term is cancelled away by a 
counterterm coming from — VPfci„ so that 



D\/w VP 

{dt+ <u>,V) <u>z= — = 55 

zt w w 
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where P ~ D^wAlnw, called osmotic pressure in Ref. [ p^ , is the net remnant of the kinetic 
pressure contribution. 

Further exploiting the kinetic lore, we can tell few words about the temperature of Brownian 
particles as opposed to the (equilibrium) temperature of the thermal bath. Namely, in view of (we 
refer to the Smoluchowski regime with t > /3~^) Pkin ~ (DP — Ij)""^ where D = we can formally 
set: 



ksTkin = [ksT - ^) < ksT . (56) 

w It 

That quantifies the degree of thermal agitation (temperature) of Brownian particles to be less than 
the thermostat temperature. Heat is continually pumped from the thermostat to the Brownian "gas", 
until asymptotically both temperatures equalize. This may be called a " thermalization" of Brownian 
particles. In the process of that "thermalization" the Brownian "gas" temperature monotonically 
grows up until the mean kinetic energy of particles and that of mean flows asymptotically approach 
the familiar kinetic relationship: J ^(< >2 ~ < u >l)dx = ksT, cf. Refs. Q for more 
extended discussion of that medium particles heat transfer issue and its possible relevance while 
associating habitual thermal equilibrium conditions with essentially non-equilibrium phenomena. 

Remark 1: Once local conservation laws were introduced, it seems instructive to comment on 
the essentially hydrodynamical features (compressible fluid/gas case) of the problem. Specifically, 
the "pressure" term VQ is quite annoying from the traditional kinetic theory perspective. That is 
apart from the fact that our local conservation laws have a conspicuous Euler form appropriate for the 
standard hydrodynamics of gases and liquids. One should become alert that in the present (Brownian) 
context they convey an entirely different message. For example, in case of normal liquids the pressure 
is exerted upon any control volume (droplet) by the surrounding fluid. We may interpret that as a 
compression of a droplet. In case of Brownian motion, we deal with a definite decompression: particles 
are driven away from areas of higher concentration (probability of occurence). Hence, typically the 
Brownian "pressure" is exerted by the droplet upon its surrounding. 

Remark 2: The derivation of a hierarchy of local conservation laws (moment equations) for 
the Kramers equation can be patterned after the standard procedure for the Boltzmann equation. 
Those laws do not form a closed system and additional specifications (like the familiar thermodynamic 
equation of state) are needed to that end. In case of the isothermal Brownian motion, when considered 
in the large friction regime (e.g. Smoluchowski diffusion approximation), it suffices to supplement the 
Fokker-Planck equation by one more conservation law only to arrive at a closed system, and 



compare with the discussion of Ref. 1 1 1 
4.2 Planar process 

Now we shall consider Brownian dynamics in the direction perpendicular to the magnetic field B , 
hence (while in terms of configuration-space variables) we address an issue of the planar dynamics. We 
are interested in the complete phase-space process, hence we need to specify the transition probability 



11 



density P {'T* ,Tt , t\~T*Q, if Oj — 0) of the Markov process conditioned by the initial data if = if o 
and = at time Iq = 0. That is equivalent to deducing the joint probability distribution 
W ^5^^^ of random vectors ~^ and previously defined to appear in the form ~^ = Tt (t) — e~^*lf o 
and i? = — f^o — Olf o respectively, cf. Eqs. (15) and (44). 

Let us stress that presently, all vectors are regarded as two-dimensional versions (the third compo- 
nent being simply disregarded) of the original random variables we have discussed so far in Sections 
2 and 3. 

Vectors ^ and it both share a Gaussian distribution with mean zero. Consequently, the joint 
distribution T4^ ^5*, i?^ is determined by the matrix of variances and covariances: C — (cy) — {{xiXj}), 
where we abbreviate four phase-space variables in a single notion of a; = {Si, S2, Ri, R2) and label 
components of x by i,j = 1, 2, 3, 4. In terms of ^ and ^ the covariance matrix C reads: 



C 



( (SiSi) (5i52> (5ii?i) {S1R2) \ 

{S2S1) {S2S2) {S2R1) {S2R2) 

{RiSi} {R1S2} {RiRi) {R1R2) 

V {R2S1) {R2S2) {R2R1) {R2R2) J 



(57) 



The joint probability distribution of 5^ and R is given by 

1 / 1 



w{s',^)=Wi^) = - (^) exp --J:c^x,x, (58) 



where c^-^denotes the component of the inverse matrix C~^. 

The probability distributions of ^ and which were established in the previous sections, deter- 
mine a number of expectation values: 

g ^ {S^Sr) = {S2S2) = ^ (1 - e-'"') (59) 
cf. Eq. (14), while {S1S2) = {S2S1) = 0. Furthermore: 

/ = = {R2R2) = 2 ^^^^2 +^2 {t + e) = 2DB{t + 6) (60) 

cf. Eqs. (20, (42), (43). In addition we have {R1R2) = {R2R1) = 0. 

As a consequence, we are left with only four non-vanishing components of the covariance matrix 
C: ci3 = C31 = ci4 = C41 = {S1R2), C23 = C32 = {S2R1), C24 = C42 = (S'2i?2) which need a 

closer examination. 

We can obtain those covariances by exploiting a dependence of the random quantities and ^ 
on the white-noise term ^ (s) whose statistical properties are known. There follows: 



Si= [ rfse"'^^*"*^ [cosa;c(*-s)^i(s)+sina;c(t-s)A2(s)] 
^0 



(61) 
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S2= [ rfse-^^*"**^ [- sin Uc {t - s) Ai (s) + cos Wc {t - s) A2 (s)] 
Jo 



^1 = / "^^ 02] 2 (1 - e-f^^*-'^ cos We - s)) + Uce-I^^'-'^ sin (i - s)] Ai (s) + 

-/3e-^(*-") sinwc (i - s) + Wc (l - e-''^*"*) cos Wc {t - s))] A2 (s) 



"^'P'+u^l 



R2 



e-^(*-") coswc 



LOre ^'^^ '^^sinw, 



{t - s)) + /Je-'^f*-*) sin We (i - s)] ^1 (s) + 
- s) + /3 (1 - e-'^(*-") coswc {t - s))] ^2 (s) . 



Multiplying together suitable components of vectors ^ and ^ and taking averages of those prod- 
ucts in conformity with the rules {Ai (s)) = and {Ai (s) Aj (s')) = 2q5ijS (s — s'), where g = ^^/?, 
i, j = 1, 2, 3, we arrive at: 



1 ft 

h = {RiSi) = {R2S2) = 2q ^ / dse-^(*-«)[/3cosWe {t - s) + 

P Jo 



(62) 



1 



and 



WcSinwc (i - -s) - /Je-'^^*"")] = q— (l - 26"'^* coswci + e'^f^*) 



k={RiS2) = -{R2Si)=2q ] ^ f dse-^'^'-'^[-psinUc{t-s)- 

P +'^c Jo 



(63) 



WcCOSUJc{t- s) -l^Jce '*)] = gr-— i 

The covariance matrix C = (cij) has thus the form 



2e-'5*sina;ei-^(l-e-2/5*) 



C = 



3 /i -fc ^ 

,9 A: /i 

h k f 

\-k h f I 



(64) 



while its inverse C ^ reads as follows: 



1 

detC 



( 


f 





-h 


k 


\ 







/ 


-k 


-h 






-h 


-fc 


9 







\ 


k 


-h 







/ 



(65) 
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where det C={fg-h^- k'^f. 

The joint probabiHty distribution of ~^ and ^ can be ultimately written in the form: 



W 



(66) 



47r2 {fg - ft2 - fc2) 



exp -- 



R 



-It +2k{y xity 



2(/5-/i2-fc2) 



(67) 



In the above, all vector entries are two-dimensional. The specific i = 3 vector product coordinate 
in the exponent is simply an abbreviation for the (ordinary i?^-vector product) procedure that involves 
merely first two components of three-dimensional vectors (the third is then arbitrary and irrelevant), 
hence effectively involves our two-dimensional ^ and 

4.3 Kramers equation and local conservation laws for the planar motion 

For the purpose of evaluating local velocity averages of the Kramers equation, we need to extract the 
marginal configuration space distribution. Let us notice that 



j W(^, lt)(f^ = w{lt) = P{-r, tiT'o, to = 0, u'o) 

where the last transition probability density entry coincides with that of Eq. (44). 
Let us introduce an auxiliary (weighted) distribution: 



(68) 



W 



w 



(69) 



27ri(/5-/i2-fc2)^''P| 2i(/g-/i2_fc2) 



where 



rrt = - {hRi - kR2, /1-R2 + kRi) 



(70) 



and we recall that = (t) — e~^*C/ (t) Tto and it = — T*o — ^ito- 

The local expectation values (compare e.g. calculations of the previous subsection) read: {ui)^ = 
J UiWdlt and (ui)^ = J ufWdlf where i = 1,2. 
By evaluating those averages we get: 



(^)7? = (("1)7? . ("2)7?) = e-^'U (t) Tto+ni 



(71) 
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and 



(«?)^ - = ("i>^ - {^^f-n = jif9-h'- k') (72) 

The Fokker-Planck-Kramers equation, appropriate for the planar dynamics in its phase space 
version, reads 

^ + itV^W = f3V^ (Wit) - [V^ X Wlt]^^^ + qV%fW (73) 

where the troublesome again (in the planar case all vectors are two-dimensional) vector product 
third component stands for 

[V^xWlf],^^ = —{Wu2)--^^{Wur) . (74) 

First two moment equations are easily derivable. Namely, the continuity (0-th moment) and the 
momentum conservation (first moment) equations come out in the form 



(75) 



and 



d d 



d 

dt[{u2)Ttw] + ^^[{<A)^w\ 
That implies 



^[("1)7? ("2)7? H 



w dri 



(77) 



w dr2 



which finally sums up to a local momentum conservation law (here, the standard vector product 
on the right-hand-side contributes its first and second components only) 



w 



kin 



(78) 
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where 'l^ kin has tensor components P^^*", and ^ kin stands for a vector whose i-th component 
is equal J2j ^^7— «ii = 1,2. Here, obviously P^'" = {{ui - (Ui)j^) {uj - {uj)^))^w and only 
diagonal entries do not vanish. Clearly Pj^'*" = {iuj)^ ~ (^i)l? 

Because of 



w. 



we can introduce again — ^-^^ — ksTkin and pass to an asymptotic regime t » = jj- Then, 
we obtain Eq. (56) to be valid in the present case as well thus quantifying an overall (magnetic field 
independent) heating process involved . 

In that asymptotic regime we have cr^ = D(3 — ^ and by employing an asymptotic form of w( /?), 
Eq. (44) we recover: 

V • P fc,„ - LD/S - - 80 

\ It J w 

together with 

K{y^)^ = -DP—. (81) 
w 

So, asymptotically {t >> j3^^) the momentum conservation law takes the form (to be compared 
with considerations of section 4.1) 

+ (ii>)^ = ^^ (82) 

However an asymptotic regime does not yet imply that the right-hand-side of Eq. (83) represents 
an acceptable "osmotic pressure" gradient contribution. We additionally need a large friction regime 
to deal with a consistent picture of a Markov diffusion process in the Smoluchowski form. Indeed, to 
reproduce a universal (Ref. [D) pressure- type functional dependence on P we must employ a suitable 
form of the diffusion coefficient. The usage of D alone to define an "osmotic pressure" implies an 
apparent failure. On the other hand, the usage of Db — D p2j^^2 as suggested by Eq. (44) leads to: 

P = I?|Alnz.^-^ = ^^ . (83) 
w 2t w 

Then however the momentum conservation law displays a supplementary scaling of the osmotic 
pressure contribution, which may trivialize (die out) only for large values of (3 (an ultimate Smolu- 
chowski regime). Namely, we have: 

a + m,v'](*),.-(^)"l£^ (84) 
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Such process is yet non-Markovian and its approximation by the Smoluchowski process becomes 
reliable when /3 is large while ujc is kept moderately small. 

Basically, the large friction regime cancels all rotational features (arising due to the Lorentz force) 
on very short time scales. If we are satisfied with the non-Markovian regime of moderate friction 
but arbitrarily varying lUc (i-e. B) then, in conformity with Eq. (78), mean flows would display 
signatures of rotation that is bound to die out asymptotically. This effect can be analyzed in by 
observing that a three-dimensional extension of the vector rft of Eqs. (70), (71) asymptotically reads 
rrt = 2^/3^(^ ~ ^o) where A comes from Eq. (3) and (T* — T*o) £ R^- In view of that, we have (cf. 
Eq. (71)) curl ^ curlrft ~ (0, 0,— |j|) and the circulation asymptotically vanishes. The effect 

can be slowed down by a suitable adjustment of uic against (3. 
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